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METHOD FOR PRODUCING A CONSTRAINT-SATISFIED 
CAM ACCELERATION PROFILE 

BACKGROUND OF THE INVENTION 
FIELD OF THE INVENTION 
[0001] The invention relates generally to methods for designing the 
profile of a cam for actuating a valve mechanism. More specifically, the invention 
relates to generation of an acceleration profile for a valve operating cam of an 
internal combustion engine, the profile satisfying a plurality of valve motion 
constraints. 

DISCUSSION OF THE PRIOR ART 
[0002] Internal combustion engines use a well-known cam shaft 
system with a plurality of cams for opening and closing various valves associated 
with individual combustion cylinders of the engine. A conventional cam-actuated 
engine valve arrangement is shown in Fig. 1. Cam 101 rotates in the direction 
shown by arrow 1 13 so as to move cam follower or tappet 103 and push rod 105 
against rocker arm 107 which, in turn, causes motion of spring biased valve 111 
in an opening or closing direction for controlling communication with cylinder 
volume 115 with an input or output conduit 113. Valve 1 1 1 is biased to a closed 
or sealed position with respect to conduit 113 by biasing valve spring 109. Zero 
degree cam angle rotation is defined as when cam nose 101a is in a vertically 
upward direction as shown in Fig. 1 wherein valve 1 1 would be in a fully open 
position. 
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[0003] At the very beginning of the cam design process, a cam 
designer may be presented with design parameters, such as overlap volume, 
intake valve closing volume, exhaust pseudo flow velocity and blow down 
volume. Additionally, manufacturing constraints such as the smallest radius of 
curvature that can be ground with a specific grinding wheel play a roll in the 
design process. 

[0004] Computerized techniques allow designers to specify how the 
valve is to move by specifying the valve acceleration. These computerized 
techniques then determine the shape the cam needs to take in order to deliver 
the desired valve acceleration profile as the cam makes a total rotation. 

[0005] Unless a design engineer is extremely lucky, the initially 
selected acceleration profile for the cam will not meet all of a plurality of valve 
motion constraints without adjusting the initial profile. Prior techniques for 
transforming draft acceleration curves into an acceleration profile that meets all 
valve motion constraints are known, wherein a plurality of scaling constants are 
sought to scale the various acceleration pulses formed by the acceleration curve 
such that the valve motion constraints will be satisfied. In known systems, there 
are four valve motion constraints but only three scaling constants due to the 
nature of the acceleration profile curve. Hence, a fourth design variable is 
chosen to be an adjustment design point acceleration value of the design 
engineer's choosing. 

[0006] The constraint satisfaction problem has conventionally been 
solved as a non-linear four-dimensional root-finding problem. The adjustment 
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acceleration value and the three scaling constants have in the past been 
adjusted by generic root-finding software in an effort to determine values of these 
four design parameters that yield an adjusted trial curve that meets all constraints 
to within an acceptable error tolerance. There are problems with this known 
approach. First, sometimes the known approach does not succeed or it does not 
deliver a highly precise solution. Secondly, this known optimization approach is 
more computationally expensive than can be tolerated during interactive design 
within many popular computing environments (e. g., Matlab / Simulink). Hence, a 
faster approach is needed. 

SUMMARY OF THE INVENTION 
[0007] In one aspect of the invention, a method for generating an 
acceleration profile for a valve operating cam of an internal combustion engine, 
wherein the profile must satisfy a plurality of valve motion constraints, begins with 
generating a valve acceleration versus cam angle draft curve by specifying a 
plurality of points of desired valve acceleration versus cam angle and using a 
curve fitting routine to form the draft acceleration curve interconnecting the 
plurality of points. A set of equations is developed, one for each of the plurality of 
constraints in terms of parameters of the draft acceleration curve and in terms of 
a plurality of scaling factors, one for each section of the draft curve between roots 
thereof. A determinant for the set of equations is formed. A point on the draft 
curve is selected as an adjustment point, and the adjustment point is varied to an 
adjustment acceleration value that forces the determinant to substantially zero. 
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The curve fitting routine is then used again to generate an adjusted acceleration 
curve which includes the adjustment acceleration value. The set of equations is 
solved for values of the scaling factors as a function of parameters of the 
adjusted acceleration curve, and sections of the draft acceleration curve between 
roots thereof are multiplied by the resultant values of corresponding scaling 
factors to generate a constraint-satisfied acceleration profile. 

BRIEF DESCRIPTION OF THE DRAWING 
[0008] The objects and features of the invention will become apparent 
from a reading of a detailed description, taken in conjunction with the drawing, in 
which: 

[0009] Figure 1 is a perspective view of a conventional cam-operated 
valve opening and closing mechanism for an internal combustion engine; 

[0010] Figure 2 is a graph of a cam acceleration profile showing an 
initial draft set of points and a continuous curve fitted among the points; 

[0011] Figure 3 is a graph of valve velocity versus cam angle resulting 
from the initial draft acceleration curve of Fig. 2 prior to adjustment of the profile 
to meet valve motion constraints; 

[0012] Figure 4 is a graph of valve lift versus cam angle resulting from 
the initial draft acceleration curve of Fig. 2 prior to adjustment to meet valve 
motion constraints; 
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[0013] Figure 5 sets forth a graph of valve velocity versus cam angle 
resulting from an acceleration curve which has been adjusted to meet valve 
motion constraints; and 

[0014] Figure 6 sets forth a graph of valve lift versus cam angle 
resulting from an acceleration curve which has been adjusted to meet valve 
motion constraints. 

DETAILED DESCRIPTION 

[0015] Suppose 1(9) defines valve lift as a function of the rotation 
angle 9 of the cam producing that lift. The second derivative of / with respect to 0 
is commonly referred to as the valve acceleration profile a(9). 

[0016] Fig. 2 shows an example valve acceleration profile for a cam, 
such as cam 101 of Fig. 1. The horizontal axis indicates cam angle. Cam angle 
zero corresponds to maximum lift — i.e., the angle where the nose of a cam lobe 
101a contacts the follower 103. Negative angles correspond to valve motion 
induced by the opening side of the cam lobe and positive angles indicate motion 
induced by the closing side of that lobe. 

[0017] The square waves 220 and 222 on the left and on the right of 
Fig. 2 are respectively called the opening and closing ramps of the acceleration 
profile. Acceleration is zero from angle -180° to the beginning of the opening 
ramp, and from the end of the closing ramp to +180°. Between the two ramps 
lies a typical valve acceleration curve, often called an acceleration profile, that is 
composed of three large acceleration pulses. These are the positive opening 
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pulse 230, the negative valve spring pulse 232, and the positive closing pulse 
234. Observe that the acceleration over the two positive pulses is always 
positive except at their boundaries, where the acceleration is zero. Similarly, the 
acceleration throughout the negative pulse is always negative except at its 
boundaries, where it is zero. For purposes of discussion throughout this 
description, it is assumed that draft acceleration curves between the square- 
wave ramps always consist of a positive pulse, followed immediately by a 
negative pulse, finally ending with a second positive pulse. There are no zero 
acceleration values except those occurring at the boundaries of the three pulses. 

[0018] In typical cam design processes, only the three pulses 230, 232 
and 234 between the two opening and closing ramps 220 and 222 are adjusted 
to create a desirable valve motion. Ramps, and their positioning within the 
acceleration profile, once set, are not typically varied. A design engineer will 
add, delete and move points that sketch out a desired acceleration curve or 
profile. A curve fitting routine, or spline, generates a curve passing through these 
points of the designer's choosing to define the cam acceleration profile a(0) 
between ramps. 

[0019] The designer's initial rough sketch 200 connects the 
acceleration data points shown as small circles in Fig. 2 such as 240, 242, 244, 
246, 214, etc. The draft acceleration profile 202 is generated by an initial 
application to the data points of a preselected spline algorithm. The data points 
are known as "knots". 
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[0020] There are four valve motion constraints that the acceleration 
profile must meet. 

[0021] The valve velocity implied by the opening ramp 220 and main 
acceleration profile must match up to the end velocity v c implied by the closing 
ramp 222-i.e., v(9 c ) = v c . 

[0022] Similarly, the valve lift implied by the opening ramp 220 and 
main acceleration profile must match up with the valve lift l c implied by the closing 
ramp 222-i.e., I(9 C ) = k- 

[0023] Additionally, the valve lift must achieve a certain maximum 
value at the nose of the cam or cam angle zero. This imposes two additional 
constraints. First, the valve lift must be some pre-selected value at cam angle 
zero (1(0) = / m ax). Secondly, the valve velocity must be zero at cam angle zero 
M0) = 0). 

[0024] As noted previously, the designer must be extremely fortunate 
to meet these constraints without adjustment of the initial draft of an acceleration 
profile. Fig. 3 is a graph of valve velocity versus cam angle where the constraints 
have not been met. Note at area 300 of the curve of Fig. 3, that the graph shows 
an end velocity of the cam which does not match up with the velocity generated 
by the closing ramp of Fig. 2. 

[0025] Similarly, Fig. 4 is a graph of valve lift versus cam angle 
resulting from an initial draft acceleration curve prior to adjustment which does 
not meet the valve motion constraints. Area 400 of the graph of Fig. 4 
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demonstrates that the valve lift generated by the draft acceleration curve of Fig. 2 
does not match up with the valve lift generated by the closing ramp of Fig. 2. 

[0026] With the acceleration profile as generally depicted in Fig. 2, the 
four constraints set forth above may be expressed in terms of parameters of the 
initial draft acceleration profile. With reference to Fig. 2, let a(9) be a draft 
continuous valve acceleration curve defined on the interval [9 0l 9J. Let 9 0 , 9 1t 9 2 
and 9 C be the only roots of a in the interval 9 0 to 9 C as shown in Fig. 2. We now 
define a new adjusted continuous acceleration function in terms of a as 



c x a(e) 9 0 <e<o x , 

c 2 a(9) 0,<0<0 2 , 

c 3 a(e) e 2 <e<%, 



[0027] c 1t c 2 and c 3 are three scaling constants to be respectively 
applied to acceleration pulses 230, 232 and 234 of Fig. 2. 



8 



Att rney R f. 706466 
HDP Ref. 0364-001622 

[0028] If a valve undergoes acceleration a(9) and has velocity v 0 and 
lift / 0 when 9 = 9 0l then the lift / c when 9 - 9 C for that valve can be shown to be 

(1) K =k- e ok +, o +c i A +c 2 I 2 +c 3 Z 3 , 
where 

[0029] Similarly, if a valve undergoes acceleration a(6) and has a 
velocity v 0 when 0 = 9 0 , then the velocity v c when 0 = 0 C for that valve is 

(2) v c =c,p;+c 2 r 2 +c3F3+v 0 , 

where 

r, = L>K 

r, 

and K 3 = \"'a{s)b. 

[0030] If a valve undergoes acceleration a(9) and, when 9 = 9 0 , that 
valve has a velocity v 0 and lift / 0 , then at 9 = 0°, that valve will have lift 

(3) l{6) = -v 0 Q 0 +l 0 +c i -L 4 + c 2 -L 5 , 

where L 4 = £ J* a{s)dsd9-e x J' 1 a(s>fc 
and I 5 = f° f a(s)dM8. 
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[0031] Finally, if a valve undergoes acceleration a(6) and, when 6 - 9 0 , 
that valve has velocity v 0 , then when 9 = 0 the valve velocity is 



(4) 



v(Q) = v 0 +c r V 1+ c 2 V 4 



where V x = f 9 ' a(s)ds, 

and V 4 = J" a(s)ds. 



[0032] It can be shown that the above four constraints can be satisfied 
if and only if the vector 6 = (d, c 2 , c 3 ) T satisfies the matrix equation 



(5) 



U L 



2 -^3 



0 



F, V A 0 



-fc-^ok-'o + O 



e oV 0 -/o+^„ 



- v ft 



[0033] Furthermore, it can be shown that a unique non-zero solution c 
to equation (5) exists if and only if 



(6) determinant 



L x L 2 L 3 



L 4 L 5 



y 3 

0 



Q n v n -L+L 



-v„ 



= 0. 



Uniqueness follows from the fact that the determinant of the lower left 3x3 
submatrix from the matrix in equation (6) above is never zero, so that the rank of 
the matrix is always 3 or larger. 
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[0034] Suppose one selects an adjustment point or knot (Oh z$ e S, 
where Q 0 <9 k < 9 C and z k * 0 (see point 244 of Fig. 2). Define the function D(z k ) 
as 



determinant 



(Afo) L 2 (z k ) Z,(z t ) -fc-0 o ]v o -/ o + / ( 
V 2 (z k ) V 3 (z k ) v c -v 0 

r.w ^(O 0 



Note that the determinant depends on a, which in turn is uniquely defined by the 
points in S that a interpolates. Thus, D can be thought of as a function of the 
non-zero interpolation value z k . For a new value of z k , D(zi<) is calculated by first 

finding the spline a that interpolates the set S , where S is the set S with the 
point (0*,z t ) replaced by (0 t ,i t ). Then entries Lj,..., L 5 and Vi,..., V 4 are 
determined from adjusted a. 

[0035] The question becomes: near z k is there a value f^for which 

D(z k ) = 0? If so, then the trial acceleration curve that interpolates the point set S 

could be replaced by the trial acceleration curve that interpolates 5. The 
resulting trial acceleration curve would look very similar to the curve that 
interpolates S (since z t is "near" z k ). It may therefore be an acceptable 

replacement for the original a. The new a will be a curve for which a scaling 
exists to solve the constraint equations developed above. 
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[0036] It should be noted that the basic goal in moving knot z k is local 

modification of the valve acceleration profile so that the determinant of equation 
(6) becomes zero. This goal may be accomplished equally well by moving two or 
more knots of the spline in concert within a localized region of the curve. 
However specifically implemented, the basic goal remains the same: add or 
subtract area from the acceleration profile locally to produce a curve for which 
equation (6) is satisfied. 

[0037] Hence, to produce a constraint satisfying acceleration profile or 
curve a from the draft curve a that meets the constraints specified above, one 
performs the following steps. 

[0038] Select a point (6 k , zk) in the set S such that z k is not equal to 

zero. 

[0039] For the function D(Zk) defined above, find a non-zero value z k 
that satisfies D(z k ) = 0. For example, one could use a root determination 

method, such as Newton's method, on the determinant. 

[0040] Replace the draft acceleration curve a with a curve generated 
by a spline using all the points of the previous curve except the adjustment point 
being replaced by (0 t ,£ t ) 

[0041] Form the matrix equation (5) and solve for the unique solutions 
to that equation for the three scaling factors c„ c 2 ,c 3 to be respectively applied to 

the acceleration pulses 230, 232 and 234 of Fig. 2. 

[0042] The new constraint-satisfied continuous acceleration function is 
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o(9)- 



c x a(e) 0 o <e<e x , 
c 2 -a(e) e x <e<e 2 , 

c 3 -a(0) 0 2 <O<Q C . 



[0043] The method discussed above assumes that a trial acceleration 
curve i(9) meets the following conditions. 

[0044] 1. a(9) is a piecewise polynomial interpolating function 
generated by the shape preserving algorithm defined below. 

[0045] 2. a(9) is a continuous valve acceleration curve defined on the 
interval [9 0 , 9 C ]. 

[0046] 3. The points 9 0 , 9 1t 9 2 and 9 C satisfy 9 0 <9 1 <0<9 2 <9 C and 
are simple roots of a. That is, these points are where the curve a is zero, and a 
is positive in the interval (9 0 , 9i), negative in (9 1t 9J, and positive in (9 2 , 9 C ). 

[0047] Below, a revised algorithm for creating shape preserving 
quadratic splines is presented. The basic algorithm is due to Schumaker, see 
Larry L. Schumaker, On Shape Preserving Quadratic Spline Interpolation, SIAM 
J. Numer. Anal., 20(4):854-864, 1983. The algorithm set forth below, like the 
unrevised version, produces continuously differentiable quadratic splines in such 
a way that the monotonicity and/or convexity of the input data is preserved. The 
revised algorithm has the additional property that the splines it produces are 
more nearly continuous in the y -coordinate values of the knots to be 
interpolated. 
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[0048] The lines of the algorithm marked with an "*" indicate where the 
algorithm has changed from the original. Input to the algorithm is a set of n 
knots (points to interpolate) { (t„z,\ i = 1, n, distinct }. 

Algorithm 1 (Schumaker-revised) 

1 . Preprocessing. 

For i = 1 step 1 until n-l, 

/ I = [fc +1 -rf+fe + ,-z i ) 2 ] W2 

5 /=( z / + i-*iM'/ + i-'/) 

2. Slope Calculations. 

For i = 2 step 1 until n-l, 

3. Left end slope. 

5,.=(38,-s 2 ) / 2 

4. Right end slope. 

^=(3S„- 1 -V,)/2 

5. Compute knots and coefficients. 

j = 0. 

For i = 1 step 1 until n-l, 

ifs i + s M =2& i 

j = j + Uj =t i ,A J = z i ,B j =s i ,C j ={s M -s ( )/2(t M +t i ) 
else 

a = (s i -S i \b = (s M -8 i ) 
if ab>0 

^={b-t M+ a-t)l{a + b) 

* elseif a = 0 

£i = _ ^ ' jiPT " (^i+l _ h ) 

m = 1; 

while ^-t M =0 

* m = 2m 

* endwhile 

* else ifb = 0 
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6-'i+C-i±rfe.i-0 



m = 1; 

while %,-t. = 0 
* m = 2m 

£> i =t i + mC > (f M -t,) 

* 



e/7c/w/7//e 
e/se/7 |a|<|6| 

e/se 

4,=^+6fc +1 -0/(-s, + .-^) 
5,. = (2* f -* w )+(* w -*,)& -0 

7 = 7 + U, = t„Aj = = s„Cj =r\,/2 

j = j + \, Xj = $,Aj = z, + * ( fe - 0+ 7/ fe, - ^,) V2, 

The following theorem can be mathematically proven and concludes that for 
every trial acceleration profile formed as a spline produced by Algorithm 1 , it is 
nearly always possible to produce a constraint-satisfied acceleration curve.. 

Theorem I. Suppose ax(t) is the shape preserving quadratic spline determined 
by Algorithm 1 for a set of knots 

where t Jy j = l,...n are distinct and increasing. When X = 0, suppose a 0 (e) is 
positive for 9 € (9 0 fi x ), negative for 9 e (0„e 2 ), and positive in 9 e {9 2 ,Q C ), 
where 9 0 <9 1 <0<9 2 <Q C . Suppose further that 
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that 0 O =/, and Q c =t n , and that for some indices i and j,t. =0, and t } =9 2 . Let 
L t ,i = 1.....5, and ^.,i = l,...,4, be defined as set forth above with a = a^. Lef 
v 0 .v c ,/ 0 ,/ c and l m be any constants such that 

-v o I 4 -^(0 o v o -/ o + / m „)^O. 
Then fhere exists at least one value of X , say X 0 , such that 



(7) lim det 



\ L 2 L 3 -[0 c -0 o ]v o -/ o +/ ( 
Vi Vi ^ v c -v 0 
L 4 L 5 0 ©oVo-Zo+U 



= 0. 



\K V, 0 -v, 

[0049] Under the hypotheses set forth in the theorem, L A , V l , v 0 , 0 O , 

/„, and l m do not depend on X. Therefore, Theorem I shows that 

whenever- v 0 Z 4 -F,(0 o v o -/ 0 +/ max )^0, the determinant in equation (7) can 

always be made arbitrarily close to zero by adjusting a properly located knot of 
the trial acceleration curve. From a computational point of view, it is nearly 
always true that only an approximate zero can ever be found to highly nonlinear 
equations, regardless of the solution technique. Theorem I in effect 
demonstrates that there is always a "numerical" solution to the constraint 
satisfaction problem. So long as -v o Z 4 -K l (0 o v o -/ o +/ max )^O, determinant (7) 

can always be made arbitrarily close to zero by adjusting X, and hence a 
constraint satisfied curve can always be produced from a trial curve that meets 
the hypotheses of Theorem I. 
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[0050] Note that while a(6) may be continuous across the roots 0f and 

02, the derivative of the constraint satisfied acceleration curve —(e) will not be. 

d0 y ' 

The derivative -^(9) is typically called the "jerk" of the valve motion. Use of the 

du 

method of this invention presupposes that a valve acceleration curve with jump 
discontinuities in the jerk at 61 and 62 is acceptable. 

[0051] Testing has been carried out on the method set forth above. So 
long as the design point (6 k , z$ (i.e., the point that is adjusted to make D(z$ = 0) 
is not too near neighboring points (9 k - 1t z k -i) and (G k+1 , z k+1 ), the following 
observations are generally true for most cases tested: 

[0052] The acceleration value z k (knot 244 of Fig. 2) need move only a 
tiny amount (see arrow 244a of Fig. 2). 

[0053] Provided 1(0) - l max is not too large, scaling constants typically 
differ from 1 by only a few percent. Therefore, the change to the trial curve is 
usually difficult to perceive. Hence, the method yields a constraint satisfied curve 
that looks quite similar to the trial curve 202. 

[0054] When the initial draft acceleration profile has been modified in 
accordance with the above method, the constraints will be satisfied as seen from 
Figs. 5 and 6. Fig. 5 shows at area 500 that the valve velocity resulting from the 
adjusted acceleration profile will match that generated by the end ramp of Fig. 2. 
Similarly, Fig. 6 shows that at area 600 the valve lift will match that required by 
the end ramp of Fig. 2. 
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[0055] To assure a solution to the nonlinear equation D(z k ) = 0 exists 
and thus assure success in meeting the valve motion constraints, the selection of 
an adjustment point should be made in accordance with the following. 

[0056] First, it is recommended that the trial or draft curve contain five 
or more distinct knots, e.g., 240, 242, 244, 246 and 214, of Fig. 2 which have 
distinct cam angle coordinates within interval [0, 0 2 ], . 

[0057] Second, the adjustment point (knot 244) should be selected 
such that the two knots immediately left (240, 242) and the two knots 
immediately to the right (246, 214) of the adjustment point 244 have cam angle 
coordinates 6 that are equal to or between zero cam angle and the third root 6 2 of 
the acceleration curve. 

[0058] These two recommendations insure that only the area of the 
design curve 202 that is between cam angle zero and cam angle 6 2 is affected by 
a change to the adjustment point z k . 

[0059] In conjunction with selecting the adjustment point in accordance 
with the above recommendations, the curve fitting routine or spline used to 
generate the adjusted acceleration profile is optimized as shown above by 
insuring that the quadratic spline will only alter the initial draft acceleration curve 
at segments between two knots on either side of the adjustment point. In other 
words, for example, if the adjustment point 244 of Fig. 2 is moved positively or 
negatively as shown by arrow 244a, the resultant adjusted acceleration profile 
generated by applying the spline to the new data set with the altered point 244 
will change the original acceleration profile curve only in segments 241, 243, 245, 
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and 247-i.e., those segments of the acceleration profile between the two points 
on either side of the adjustment point. 

[0060] The invention has been described in connection with an 
exemplary embodiment and the scope and spirit of the invention are to be 
determined from an appropriate interpretation of the appended claims. 
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